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1. Introduction 


Mass Reweighting [[!]] is an interesting and efficient method to correct and to include effects of 
quark masses. It can be used for tuning, e.g. the strange quark mass m s in a 2+1 simulation or the 
isospin splitting in the up-quark mass m u and the down-quark mass m d . Moreover it can be applied 
to understand the mass behavior of observables. For small corrections it is applicable and more 
efficient than new simulations. Mass reweighting involves the evaluation of fermion determinants 
which can be rewritten by an integral representation. This integral representation can be estimated 
by Monte Carlo integration which needs around 100 inversions of the Dirac operator to control the 
stochastic noise efficiently. 

The reweighting factor enters the measurement of an observable by [|2|] 

(ov = ^ = {ow) (1.1) 

where the mass reweighting factor for nj flavors of quarks 


» , =n 


■' det D„ 
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is normalized with W = W/(W). Here, the Dirac operator is given by the 
Dirac operator D m = D + m. The reweighting factor can be rewritten as 
matrix M with 

W = -=-. 

det rijii det M 


clover improved Wilson 
a determinant of a ratio 

(1.3) 


In general lattice simulations are done in the isospin symmetric limit in the light quark sector by 
setting the light quark masses to the average light quark mass m ud = ().5(m u + ni d ). The idea is to 
use mass reweighting to introduce isospin breaking. The reweighting is performed by splitting up 
the light quark masses by keeping the average quark mass constant 2 m ud = m u + m d = const and it 
follows with the mass shift Any = m d — m u 


m u = m ud - 0.5 • Am ud +- m ud 
This leads to the isospin reweighting factor 

W„ 0 = 
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Now, the determinant of the non-hermitian matrix M can be rewritten by an integral representation 
given by 

^^ = /D[T]]exp{-T7 + MT7} (1.6) 

which holds for X(M + M') >0 [^]] and the normalized integral measure is given by D[ 77 ] = 
n”_i dxjdyj/n with = Xj + iyj. The integral eq. ( |1.6[ ) can be estimated stochastically 
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Table 1: The table shows the analyzed CLS - ensembles generate with rif = 2 - ff{a) improved Wilson 
fermions with m ur ] = m u = |Q|. We used ensembles of three different lattice spacings a with pion masses 


m K from 580 MeV down to 192 MeV and lattice volume V/a 4 . The number of configurations N cn f g are 
seperated by MDUs/config with R act the relative number of active links. The maximal reweighting range is 
given by the quark mass shift in the /V/.S'-scheme A niR max . The renormalized mass is defined as in [ 0 ]. 

ID 

V/a 4 

a [fm] 

m n [MeV] 

Ncnfg 

MDUs/config • R act 

Nm R , max [MeV] 

A5 

64 x 32 3 

0.076 

330 

202 

20-1 

4.43(60) 

E4 

64 x 32 3 

0.066 

580 

100 

16-0.37 

7.1(17) 

D5 

48 x 24 3 

“ 

440 

503 

8-1 

5.9(10) 

E5 

64 x 32 3 

“ 

440 

99 

160-0.37 

6.01(96) 

F7 

96 x 64 3 

“ 

270 

350 

16-0.37 

5.01(38) 

G8 

128 x 64 3 

“ 

192 

90 

8-1 

5.86(56) 

07 

96 x 48 3 

0.049 

270 

98 

40-1 

5.63(40) 


by drawing Nj, pseudofermion fields r) distributed via the normalized function °< exp{— r/ ' r/j and 
I the unit matrix with dimension of M. Note for every drawn field rj n y- inversions of the Wilson 
Dirac operator D m have to be performed. 

In general mass reweighting introduces fluctuations which increase the statistical error. These 
fluctuations are the ensemble fluctuations, introduced by the ensemble average in eq. ( PI ), and 
the stochastic fluctuations, introduced by the stochastic estimation of the integral eq. The 

fluctuations are given by the variance averaged over the ensemble and the pseudofermions rj. We 
will define the variance of the integral representation eq. ( p| ) by a 2 = ((vvvv 2 )^} — ((vv) T) )((vv t )^) 
with the stochastic estimate w(U, tj) = 1 /V, ; £ ; exp{— rfj ( M—l)r\i] . By performing the rj-average 
()jj (i.e. all T], independently) the fluctuations are given for finite /V () by 


Nr, (det(M + Mt-/)) + Nr, (detMM+) ( detM ) ( detAft) 


( 1 . 8 ) 


which holds for X(M + M 1 —I) >0. The stochastic fluctuations for one configuration are given by 
neglecting the ensemble average () and vanish for Nj, —>• °o. Moreover by introducing a mass inter¬ 
polation between the start and the target mass the stochastic fluctuations can be further controlled, 
i.e. if no Dirac operator has a zero eigenvalue during the interpolation the condition ) > 1 

can be insured [|5]]. In this case the number of inversions is N inv c< N -Nr,, where N is the number of 
interpolation steps. Note for many reweighting cases it is effecient to use the even-odd precondi¬ 
tioned Wilson-Dirac operator (e.g. see [[(]]), however we do not find an improvement in the case of 
isospin reweighting. The ensemble fluctuations can be tamed by including additional quarks into 
the reweighting process, e.g. in the case of isospin reweighting the fluctuations are minimized by 
keeping the average quark mass m uc i = 0.5(m u +m l j) constant during the reweighting. 

Here, we will discuss mass reweighting by introducing an isospin breaking in the light quarks. 
We will show the scaling of the different fluctuations (ensemble and stochastic) and how the up- 
and down-quark mass can be extracted from the analyzed ensembles (see tab. [I]). 
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Figure 1: The relative stochastic fluctuations of the isospin reweighting factor is shown as a function of the 
renormalized quark mass, o 2 r is the average of the stochastic variance estimated using Nq = 6. 


2. Isospin Reweighting 


The isospin reweighting factor eq. ( |1.5| ) can be expanded in A rrr ud 

Wise = dCt q e m t 'of Pm " = 1 + Am ud Tr ( D 7n ud ) + ®( Am \d) (2- 1 ) 

m ud 

by using detM = exp{Tr(ln(M))}. The fluctuations a 2 in eq. ( |T 8 | ) of the isospin reweighting 
factor can be expanded in A m 2 ud . It can be shown that the stochastic fluctuations o 2 t decouple 
from the ensemble fluctuations c 2 m with o 2 = G 2 + o 2 ns . The stochastic fluctuations in the isospin 
reweighting case are given by 


<*l (Njnv) \ 

w 2 / 


At A / Tr_1__ 

N,m \ 




( 2 . 2 ) 


The ensemble fluctuations are 


(W ) 2 


= Am td var 




(2.3) 


with the variance var(O) = (O 2 ) — (O') 2 . Note this is true because the Dirac operator is 75 - 
hermitian. Now, the cost can be derived by demanding that the stochastic fluctuations do not 
dominate the ensemble fluctuations, e.g. that c 2 t (Ni nv )/ (j- ^ ^ 0 . 1 . 

By using the analyzed ensembles, listed in table [I], it is possible to deduce numerically the 
scaling behavior of the fluctuations in the quark mass and the volume. However in the case of 
the stochastic fluctuations the trace of the Wilson Dirac operator is known in chiral perturbation 
theory, e.g. as in 0, by (Tr CD£Ft> “ if’ The numerical analysis is consistent with this behavior. 
It follows for the stochastic fluctuations (see fig. [j]) 


Am 4 R V 1 

M. m J y 
iy inv rn R 


(2.4) 
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Figure 2: The scaling of the ensemble fluctuations for different volume behaviors is shown. The left figure 
shows the quark mass behavior for a volume behavior of \/V and the right figure for a volume behavior of 

V 3 / 4 . 


by using the scale ro 0] to form dimensionless quantities. By fixing the volume behavior to V 
we perform a fit (black, solid) to the quark mass behavior by using the ensembles E4 (green, 
triangle), D5 (red, diamond), E5 (black, triangle), A5 (cyan, circle), 07 (magenta, diamond), F7 
(blue, square) and G8 (green, star). For each ensemble we compute o; t and oy ns for two values of 
Am = Am max /2.Am max (see tab. |]). The quark mass behavior is given by r' = 2.63(5). For the red 
lines the quark mass behavior and the volume behavior is fixed to V/m q and only ensembles with 
pion masses < 340MeV are included. The data show a good agreement with the expectation from 
chiral perturbation theory for pion masses < 340MeV. 

The leading term of the ensemble fluctuations eq. (2.3) is proportional to var(TrZ) 2 ). Numer¬ 
ically we observe a weak volume dependence V q with q < 1. Similar to the stochastic fluctuations 
the ensemble fluctuations can be written as 

2 Am 4 R V q 


1 


m'r 


(2.5) 


In general the simultaneous deduction of the volume and the quark mass behavior is difficult. 
A varied volume behavior changes simultaneously the mass behavior. A good fit is given for a 
volume scaling of q = 0.25 (see left figure |j) which gives a mass behavior of r = 3.85(13) for all 
ensembles (black line) and r = 3.94(31) for ensembles with pion masses < 340MeV (red dashed 
line). A weaker volume behavior is also supported by comparison of D5 and E5 ensembles, which 
gives q ~ 0.46. However by assuming a similar quark mass behavior as in the case for the stochastic 
fluctuations with r = 3 the scaling in the volume is roughly given by q ~ 0.75. In the right figure 
|2] we fixed the volume behavior to q = 0.75 which gives a mass behavior of r = 2.83(13) by 
including every ensemble (black line) and r = 3.04(31) by including ensembles with pion masses 
smaller than < 340MeV (red, dashed line). We conclude that the volume behavior is given by 
q ~ 0.25 ... 0.75 by a simultaneous variation of the quark mass behavior from r ~ 4... 3. 

The cost of isospin reweighting can be estimated from the ratio 


o, t (N illv )/ a, 


2 Kt ( Lm n) 

- 

ens 


k'ens Ninv ' T) 


with 


V 
K st 
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= le-3 


(2.6) 
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Figure 3: The figures show the average quark mass m ll( ] (left) and the mass splitting mj — m u plotted versus 
Rj, after fixing R\ and Ri. The masses are renormalized in the MS-scheme at 2GeV. 


for q = 0.25 and r = 4. For the G8 ensemble follows N- inv ~ 200 for a ratio of 0.1. 


3. Quark Masses 


The continuum limit can be performed on a line of constant physics. This line can be defined 
by keeping dimensionless ratios of physical quantities constant. These fix the bare mass parameters, 
here, a quenched strange quark with m s , the isospin mass splitting A m uc i and the average light quark 
mass mud- We take the ratios 


0.5(m|o +m|±) 
(0.5(/ x o+/ JC± )) 2 


r 2 


o - m\ ± 

(0-5(//f 0 + f K ±)) 2 


and /?3 


(0.5(f K o+fK±)) 2 


(3.1) 


with the meson masses, the pion m K ±, the neutral kaon m K o and the charged kaon m K \ and the kaon 
decay constants f K » and f K ±. The physical values of the ratios are taken from [ 0 1 and we assume 
0.5 (f K o + [k i ) = 155 MeV. Now, the strategy is to use R\ to fix m s , which is done in [j4| and Ri 
to fix the isospin splitting Amy. Afterwards the ratio A? 3 is used to extrapolate the light quarks 
towards the physical limit. 

We measure the PCAC mass on the analyzed ensembles (see tab. [j]) and convert them into 
the MS-renormalization scheme. The dimensionless ratios IF and R 2 are given in the lowest order 
chiral perturbation theory up to 0(Am 2 d ,m^ d ) by IF = jj Amy( I + Cmy) and R 2 = jj m uc i with 
constants B , C and F 2 . Now, it is possible to perform extrapolations towards the physical point in 
the light quark masses. For the average light quark mass this is shown in the left figure of fig. || by 
assuming m uc i(RT,) ~ a\R 2 . By using the F7 and G 8 ensemble the average light quark mass at the 
physical point at finite lattice spacing of a = 0.066 fm is given by m uc ip = 3. 19(11) MeV. For the 
mass splitting in the light quarks (see right plot in figure |3]) we assumed Am ut: j(R 2 ) ~ bo + b\R 2 . By 
using the data of the E5, F7 and G 8 ensemble it follows for the mass splitting Any = 2.49(10) MeV 
at finite lattice spacing a = 0.066 fm. 

The isospin effects enter the observable by the isospin reweighting factor which scales propor¬ 
tional to Am 2 d . In the case of the PCAC mass the statistical error is too big compared to the effect 
of the isospin reweighting correction. A determination of this effect is only possible for larger 


6 














Isospin Effects by Mass Reweigliting 


Jacob Finkenrath 


statistics. Neclecting the sea quark effects proportional to A m~ td , by setting the isospin reweighting 
factor to one, the isospin quark mass splitting is given by A m u d = 2.52(10) MeV. However the 
isospin reweighting effects increases for smaller quark masses and we want to reduce the statistical 
error to figure out the isospin effects for example in the pion mass. In order to perform a continuum 
limit the statistics has to be increased and other ensembles have to be included. 


4. Conclusion 


Isospin mass reweighting needs a moderate numerical effort. The analysis shows that the cost 
scales with ( LMps) 2 L for a volume scaling of the ensemble fluctuations with \/V and is around 
200 inversions of the Dirac operator for the G8 ensemble which has a pion mass of 192 MeV at 
a volume of V/a 4 = 128 x 64 3 . By using the introduced dimensionless ratios R\, Rj and R$ it is 
possible to extract the light quark masses. The isospin mass spliting is A m U( j = 2.49(10) MeV and 
the average quark mass is m uc i,R = 3.19(11) MeV at finite lattice spacing of a = 0.066 fm. Although 
a more careful analysis is needed to extract competitive numbers it shows that the tuning conditions 
are suitable to extract the light quark masses. In order to extract continuum physics the statistics 
has to be improved and QED-effects have to be included. A software package for mass reweighting 
[ 10] (see also JTT|]) is publicly available in the framework of the openQCD code[[l2|]. 
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